Strong Coupling of a Spin Ensemble to a Superconducting Resonator 
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We report the realization of a quantum circuit in which an ensemble of electronic spins is coupled 
to a frequency tunable superconducting resonator. The spins are nitrogen- vacancy centers in a 
diamond crystal. The achievement of strong coupling is manifested by the appearance of a vacuum 
Rabi splitting in the transmission spectrum of the resonator when its frequency is tuned through 
the nitrogen- vacancy center electron spin resonance. 



Building a quantum computer requires engineering a 
system that can reliably store quantum information and 
process it through a succession of quantum gates. Quan- 
tum bit implementations based on individual microscopic 
systems such as atoms, photons, electron or nuclear spins, 
benefit from a natural decoupling from environmental 
noise, which results in long coherence times [TJ; on the 
other hand, superconducting qubits |2J are macroscopic 
artificial atoms that couple strongly to electromagnetic 
fields, allowing faster single- and two-qubit gates, al- 
though with shorter coherence times [3]. It is thus ap- 
pealing to take the best of both worlds by combining arti- 
ficial and natural quantum systems in "hybrid" quantum 
circuits that would exhibit long coherence times while al- 
lowing rapid quantum state manipulation. Whereas the 
coupling strength g of an individual microscopic emitter 
to one electromagnetic mode is usually too weak to al- 
low for coherent exchange of quantum information, the 
coupling strength of an ensemble of N such systems is 
enhanced by \ /f N @HB], allowing to reach the strong 
coupling regime g\^N ^> ft, 7, where n and 7 are the 
resonator and emitter damping rates. In that perspec- 
tive, it has been proposed to couple an ensemble of cold 
molecules [7J, atoms [8 , 9j, or electron spins pTU| [TT] to 
a superconducting resonator that would mediate their 
interactions with one or a few superconducting qubits. 
This system would work like a genuine quantum Turing 
machine, with the ensemble providing a quantum mem- 
ory [T2HT3] with long storage time, and the qubits pro- 
viding the "hardware" to perform quantum gates [15J. 

Among all microscopic systems that can be coupled 
to superconducting circuits, negatively charged nitrogen- 
vacancy centers (N- V) in diamond are particularly at- 
tractive |16| . They consist of a substitutional nitrogen 
atom and an adjacent vacancy having trapped an addi- 
tional electron; their electronic ground state has a spin 
S — 1 with the states m$ = and m$ = ±1 separated by 
~ 2.87 GHz in zero magnetic field [17]. The coherence 
time corresponding to this transition has been demon- 
strated to be as long as 2 ms at room-temperature for 
samples isotopically enriched in 12 C [18J. Compared to 



atoms, N- V centers are perfectly compatible with super- 
conducting circuits, because they do not require chal- 
lenging trapping techniques or large magnetic fields to 
bring them in resonance at GHz frequency with the cir- 
cuit. Finally, in addition to the electron spin resonance 
(ESR), N- V centers have two very interesting internal de- 
grees of freedom: their narrow optical resonance might 
be used for coherently converting microwave into optical 
quantum states of the field [19J, and their coupling to the 
nitrogen atom nuclear spin could give access to coherence 
times much longer than with electron spins [20J. 

In this letter we report the observation of a vacuum 
Rabi splitting [4 , 6 , 21] [22] in the transmission spectrum 
of a superconducting resonator magnetically coupled to 
an ensemble of ~ 10 12 N- V centers with a collective cou- 
pling constant g ens /2ir = 11 MHz. This demonstrates 
the strong coupling of a microscopic spin ensemble to a 
macroscopic circuit. In a parallel experiment, Schuster 
et al. [23] observe a similar coupling to electronic spins 
in ruby and diamond (see also [!!])• The sketch of our 
experiment is shown in Fig. [T] A diamond single-crystal 
(3 x 3 x 0.5 mm 3 ) containing the N-F centers is glued 
on top of a half- wavelength niobium coplanar waveguide 
resonator, with a distance to the silicon substrate less 
than ~ 0.5/im to ensure a maximum spin-resonator cou- 
pling. The diamond is positioned in the middle of the 
resonator where the magnetic field is maximum, with its 
(001) crystallographic plane facing the chip. The spin 
Zeeman splitting can be tuned with a magnetic field B^v 
parallel to the sample surface along the [100] axis within a 
few degrees. An array of four superconducting quantum 
interference devices (SQUID) is inserted in the resonator 
central conductor, away from the diamond crystal, to 
make its frequency uj r (&) tunable with the magnetic flux 
<E> threading the SQUID loops [25 . This flux is gener- 
ated by passing current through an on-chip wire so that 
the resonator can be brought in resonance with the spins 
without changing their Zeeman splitting. The resonator 
transmission IS21I (cj) is measured with a network ana- 
lyzer, at powers low enough for the current through the 
SQUID to stay well below its critical current, so that the 
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resonator behaves linearly [25J. 
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Figure 1: Scheme of the experiment, (a-c) Photograph, top 
and side view of the sample : a diamond crystal is glued 
on top and in the middle of a superconducting coplanar res- 
onator with its (001) surface facing the chip. A SQUID array 
(d) is inserted in the central conductor to tune the resonator 
frequency cj r (<E>) with an on-chip wire producing a flux <3> in 
the SQUIDs. A magnetic field Bnv is applied parallel to the 
[100] direction, (e) Sketch of a N- V center with its vacancy 
(V) and nitrogen atom (N), as well as three neighboring car- 
bon atoms. Four equivalent < 111 > directions exist for the 
N-V axis, all making the same angle 6 = 55° with [100], and 
thus with Bnv- (f) Energy diagram of the N-V center, with 
the \ms = 0) ground state separated by uj± from the excited 
states |d=), which are linear combinations of the pure spin 
states \ms — +1) and \ms = — 1). 

We first characterize separately the coplanar resonator 
and the diamond crystal. A typical spectrum of the 
resonator alone, cooled at 40 mK, is shown in Fig. [2^l. 
The periodic dependence of the resonator frequency with 

is shown in Fig. |2j), in excellent agreement with the 
model described in [25] for parameters close to design val- 
ues. On top of the modulation curve, the quality factor 
Q = 5 • 10 4 results from the combination of the coupling 
to the external line (Qext = 10 5 ) and from internal losses 
(Qint = 10 5 ). As already observed in similar samples, Q 
decreases when the resonator is tuned to lower frequen- 
cies, possibly due to flux noise or losses in the Josephson 
junctions [25 J. Around the N- V centers frequency of 2.87 
GHz, we find Q ~ 2 • 10 4 corresponding to a resonator 
energy damping rate k = co r /Q ~ 0.9 MHz. 

The diamond crystal (of the high pressure high tem- 
perature lb type, with a nominal 100 ppm nitrogen con- 
centration) was irradiated with 2.5 MeV protons at a 
dose of 5 • 10 16 cm -2 in order to create vacancies, and 
subsequently annealed at 900 °C for 10 h to form neg- 
atively charged N- V centers. The resulting N- V dis- 
tribution is homogeneous over a 30 jam depth from 
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Figure 2: (a) Resonance line of the coplanar resonator with 
no diamond on top, measured at 40 mK. (b) Measured (dots) 
and fitted (dashed line) central frequency uo r as a function of 
the applied flux <E>, yielding a SQUID critical current 2/co = 
2.2 jiA and a bare resonator frequency cj r o/27r = 3.012 GHz. 
(c) Optically detected ESR of the N- V centers in the sam- 
ple, in zero magnetic field and at room-temperature. The uj± 
lines are fitted by a sum of two 6 MHz wide Lorentzians sep- 
arated by 10.4 MHz. (d) Measured (dots) and fitted (dot- 
ted lines) dependence of uj± on Bnv, yielding a zero-field 
splitting D/2tt = 2.873 GHz and a strain-induced splitting 
E/2tv = 4.3 MHz. 



the irradiated surface and thus over the spatial exten- 
sion of the resonator magnetic field (characteristic de- 
cay length of about 10 jam). The N-V concentration 
p = (1.2 ± 0.3) x 10 6 /im" 3 was measured by compar- 
ing the sample photoluminescence under a laser excita- 
tion at 532 nm to the photoluminescence of an individ- 
ual center in the same conditions. A typical optically 
detected ESR spectrum is shown in Fig. ^p: two reso- 
nance lines are observed even in zero magnetic field, due 
to residual strain in the crystal that lifts the degener- 
acy of states \ms = +1) and \ms = — 1), yielding new 
eigenstates |+) and | — ). The ESR frequencies measured 
as a function of £?nv [ see Fig. (2jd)] are in quantitative 
agreement with the expected transitions uj±(Bnv) be- 
tween the energy eigenstates of the spin Hamiltonian 
H NV = hDS 2 z + HE(Sl - S\) + £ N v/i B S • B NV El, 
with a zero-field splitting D/2ir = 2.873 GHz, a strain- 
induced splitting E/2tt = 4.3 MHz, a N-V Lande factor 
#nv = 2, and \i^jh = 14 MHz/mT. Here, Z denotes 
the N-V axis [see Fig. [lje)], while X and Y are the di- 
rections of nonaxial strain in the diamond matrix. Note 
that only two lines are observed in our geometry because 
the four possible N- V crystalline orientations are all at 
the same angle = 55° with the [100] direction and thus 
with Bnv- The measured ESR lines have &j/tt ~ 6 MHz 
FWHM linewidth, compatible with the expected broad- 
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Figure 3: Transmission \S2i\ of the resonator with diamond on top (at 40 mK) as a function of $ for 5nv = 0, 0.99, and 
1.98 mT. Two anticrossings are observed symmetrically around 2.88 GHz, when the resonator frequency is resonant with the 
N-V transitions |0) — ► |+) and |0) —> \— ). Red solid (yellow dashed) lines are fits to the eigenfrequencies of the coupled 
(uncoupled) resonator-spins system as described in the text. A transmission spectrum (white overlay) is also shown in linear 
units in the middle of one of the two anticrossings, for each Bnv value. Right panel: Experimental frequencies Cb± (dots) 
resulting from the fits of IS21I (cj, <£) for several Bnv values. The data points are then fitted (solid lines) to the eigenfrequencies 
of Hamiltonian H^v yielding D/2tt = 2.878 GHz and E/2tt = 7.2 MHz. For comparison, the ESR frequencies w±(5nv) 
measured at room-temperature (see Fig. |2l) are also shown as dashed lines. 



ening due to dipolar interactions with the neighboring 
nitrogen electronic spins (S = 1/2), given the nominal 
100 ppm nitrogen concentration [27J. 

We now present measurements of the resonator trans- 
mission with the diamond crystal on top. The microwave 
power used corresponds to a maximum intracavity en- 
ergy of about 85 photons at resonance and is well below 
the value required to saturate the spins. Whereas opti- 
cal pumping is required to observe the ESR at room- 
temperature, cooling down the sample to 40 mK al- 
ready provides a strong spin polarization in |ras = 0). 
Two-dimensional plots of the transmission spectrum as 
a function of are presented in Fig. [3] for several val- 
ues of £?nv, showing two avoided crossings when the res- 
onator is tuned through the N-V center ESR frequen- 
cies. To account for these observations, we model the 
resonator-spins system |4J by the Hamiltonian H/h = 
u r ( ) a t a + u) + b\_ b+ + 00 _ b L b _ + ( g+ a t b+ + g _ a t b _ + h . c) , 
whose justification is given below. It describes the cou- 
pling of the resonator mode (annihilation and creation 
operators a and a"*") to two spin- wave modes (annihila- 
tion and creation operators b± and b±) representing the 
two ESR resonances at frequencies uj±(Bnv), with cou- 
pling constants g±. The transmission spectrum |52i|(^) 
is fitted for each $ by a sum of Lorentzian peaks, whose 
central frequencies are then fitted to the eigenfrequen- 
cies of H (red solid lines in Fig. [3|. At zero field, this 
fit yields #+/2tt = #_/2tt = 11.0 ± 0.5 MHz and a co- 
operativity parameter C = # 2 /(^7) — 27, which brings 
further confirmation that our experiment is in the strong 
coupling regime. The Cj±(Bnv) data points obtained and 



shown in Fig. [3] agree very well with the eigenfrequencies 
of the spin Hamiltonian H^y for E/2ir = 7.2 MHz and 
D/2tt = 2.878 GHz. This establishes that the anticross- 
ings observed are indeed due to the magnetic coupling 
between the coplanar resonator and the N-V centers. 
The small changes in E and D with respect to the room- 
temperature values are likely due to thermal contraction 
upon cooling [28] . 

The coupled oscillator model also predicts that the 
widths of the two peaks at each anticrossing should be 
equal to the average of the resonator and spin linewidth. 
While this is approximately true at B^v — 1-98 m T, 
where the two lines have comparable linewidths of ~ 
7 MHz in the middle of the anticrossing, this is clearly 
not the case at £?nv = mT, where the line inbetween 
the two anticrossings is much narrower than the other 
lines and narrower than the coupled model predicts (see 
Fig. [3]). A more refined model is needed to understand 
this phenomenon. We also note that the 1.5 MHz FWHM 
resonator linewidth (corresponding to k = 9.4 MHz) in 
the presence of the diamond far away from the anticross- 
ings is much larger than in Fig. |2^i, indicating that this 
particular diamond sample introduces losses of unknown 
origin. 

We now justify our model for the interaction of the 
spin ensemble with the resonator. A single N- V center 
with spin operator Sk located at position rk is coupled to 
the resonator magnetic field B r (rk) by the Hamiltonian 
#k = #NvMBB r (r k ) • Sk- Here B r (r k ) = 5B (r k )(a + a 1 ") 
with £B (rk) the rms vacuum fluctuations of the mag- 
netic field at rk. Restricting ourselves to one of the two 
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|0) —> |=b) ESR transitions, this Hamiltonian can be put 
under a Jaynes-Cummings form = #k(rk)(&cr+,k + 
fl^-,k), with # k (r k ) = (^NvMB/^)^B (r k ) • (0| S k |±) 
and <j +) k ( cr -,k) the raising (lowering) operator of the 
corresponding transition [29J. The coupling Hamilto- 
nian of the ensemble of N spins to the resonator is 
H ens = a gk( r k)o"+,k + fr-c, indicating that one ex- 
citation in the resonator mode is coupled to a well- 
defined coherent superposition of spin excitations (a 
spin- wave) [TT| US] . In the limit where the number of 
system excitations is small compared to A/", this spin- 
wave behaves as a harmonic oscillator described by 
an annihilation operator b = XX#k(rk)/#ens)0"-,k> and 
the Hamiltonian can be rewritten H ens = g ens (ab^ + 
a)b) with g ens = (J pdr \g(r) J 2 ) 1 / 2 . For a homo- 
geneous distribution p and a resonator and crystal 
of lengths L and Z, this collective coupling cons tant 
can be expressed as g ens = #nv \/r]apoftw Y (&) p/2h. 

Here, rj = (1/L) ^-l)/2 ^ T dx and a = 
f |£B (r)| 2 sin 2 cp(r)dr/ f |(5B (r) | 2 dr are dimensionless 
factors describing respectively the fraction of the res- 
onator mode volume occupied by the spins and their aver- 
age orientation with respect to the resonator microwave 
field, (p(r) being the angle between the N-V axis and 
B r (r). We stress that apart from 77 and a, the ensem- 
ble coupling constant only depends on the spins density. 
For our sample, we estimate 77 = 0.29 and a = 0.81, 
which combined with the measured N-V density yields 
g ens /27r = 11.6 MHz, in agreement with the fitted value. 
This confirms that the spin ensemble is highly polarized 
and thus not far from thermal equilibrium. 

Although the strong coupling regime was reached in 
this sample, the resonator and spin linewidths should 
be reduced by 1 order of magnitude in order to imple- 
ment a quantum memory. This will require eliminating 
the extra microwave losses caused by the diamond crys- 



tal which are presently limiting the resonator linewidth. 
The spin resonance linewidth on the other hand is be- 
lieved to be limited by dipolar interactions with neigh- 
boring nitrogen electronic spins; improving the nitrogen 
to N- V conversion rate would reduce the N- V linewidth 
while maintaining the large collective coupling constant 
measured here. The very same setup would then allow 
one to investigate the coherent oscillations between the 
spin ensemble and the resonator, using a rapid tuning of 
the resonator frequency [25] in and out of resonance with 
the spins. This would open the way to the storage and 
retrieval of a given microwave field at the single-photon 
level in the spin ensemble. 

In conclusion we have observed vacuum Rabi splittings 
in the transmission spectrum of a superconducting copla- 
nar resonator magnetically coupled to an ensemble of N- 
V centers, with a collective coupling constant as large 
as 11 MHz. The position of these anticrossings is in 
excellent agreement with the electron spin resonances, 
and the measured coupling constant is well understood. 
These results therefore constitute an experimental evi- 
dence for the coherent coupling of a spin ensemble to a 
superconducting circuit, an essential step towards the im- 
plementation of more complex hybrid quantum circuits 
in which superconducting qubits, electron and nuclear 
spins, microwave and optical photons would be coher- 
ently coupled. 
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